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HI holes in galactic disks: Tracing the dark matter distribution 
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Abstract. Multiple SN explosions in disk galaxies efficiently evacuate gas and form cavities with the sizes and 
shapes of the surrounding envelopes determined by the total amount of injected energy and by the initial gas 
distribution. Such cavities are seen as HI holes when observed in face-on galaxies. Gas hydrodynamics simulations 
are performed to obtain the quantitative characteristics of HI holes that could serve for the determination of the gas 
vertical scale height and the corresponding dark matter content and its distribution. Among these characteristics 
is the ratio of the maximum column density in the HI ring surrounding the hole to the background HI column 
density and the vertical expansion velocity of gas in the HI ring. We show that in some cases the extragalactic 
background ionizing radiation may produce HI holes in the outer regions of galaxies, and can account for the 
existence of HI holes in nearby face-on galaxies with the apparent lack of an underlying stellar population. 
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I. Introduction 

The interstellar medium (ISM) shaped by the energy re- 
lease from exploding massive stars provides the best sam- 
ple for studying star formation processes on a "micro- 
scopic" level. The diversity of HI structures and associated 
Ha distributions discovered in our own and nearby galax- 
ies shows that star formation processes and dynamics of 
interaction of massive stars with the ISM are strongly en- 
vironmentally dependent. Prominent examples are found 
in the Milky Way (Heiles 1984), M31 fBrinks lTrMl Brinks 
& Bajaja[THll, M 33 (Deul & den Hartog lT330|) . Holmberg 
II (Puche et al. 11335)1 . NGC 891 (Dettmar & Schulz lT332|) . 
M101 and NGC 6946 rKamnhuis lT333l> . NGC 253 (Sofue, 
Wakamatsu & Malin [T33Ijl . SMC (Stavely-Smith ct al. 
IT337|l and LMC (Kim et al. lT333|> . IC 10 (Wilcots & Miller 
IT3381, IC 2 574 (Walter & Brinks 1133311 . Holmberg I (Ott 
et aLHED, and others. Edge-on (such as NGC 891, NGC 
4631, NGC 253) and face-on (such as Holmberg I and 

II, IC 2574, SMC and LMC) galaxies show these "mi- 
croscopic" structures in vertical and horizontal directions, 
respectively, and develop our understanding of not only 
the interrelation between the observed HI/Ha structures 
and star formation activity (Dettmar [1992 Rand 1996 
Sofue, Wakamatsu & Malin lT331 Rossa & Dettmar 2003, 
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Stewart et al. 1211331 Walter & Brinks 1999 ), but also the 
factors governing the star formation process itself. 

Among these factors, the dark matter distribution in 
radial and vertical directions is of principle importance. 
At the same time, the influence of the dark matter dis- 
tribution on the dynamics of HI and Ha structures in 
edge-on and face-on galaxies was left aside general dis- 
cussion. Only recently Ott et al. (2001) paid closer atten- 
tion to the relevance of HI dynamics in the dwarf irreg- 
ular galaxy Holmberg I to the vertical gas distribution. 
However, available estimates of the vertical scale heights 
are still based on qualitative arguments and assumptions 
on the contribution of the dark matter to the total grav- 
itational potential. In this paper we address the question 
of whether the observations of morphologies and quantita- 
tive characteristics of the HI holes in face-on galaxies can 
provide us with firm conclusions about the vertical scale 
height of the HI distribution. For this purpose, we study 
numerically the characteristics of HI flows associated with 
multiple SN explosions. We determine the dependence of 
these characteristics on the vertical scale height of the HI 
distribution and, consequently, on the vertical shape of 
the gravitational potential and dark matter distribution. 

The paper is organized as follows. In Sect.|21the numer- 
ical model is formulated, and the quantitative character- 
istics of HI holes are numerically investigated for different 
vertical scale heights of gas distribution in face-on galax- 



2 



E. I. Vorobyov and Yu. A. Shchekinov: HI holes in galactic disks 



Table 1. Model parameters 





PdO 


ro 


p* 


h, 


«* 


M d 


0.5 


0.07 


0.227 


0.02 


0.3 


7 


8 x 10 7 


0.37 


0.015 


2.26 


0.02 


0.3 


7 


3.7 x 10 s 


0.21 


0.13 


1.5 


0.03 


0.3 


8.5 


2.2 x 10 9 



^ all scales are in kpc, densities in Mq pc -3 , stellar velocity 
dispersion (v f ) in km s , the mass of the dark matter 
(Md) in Mq is computed within the radius of 2 kpc. 



ies. In Sect.[3]the influence of the extragalactic background 
ionizing radiation on the vertical shape of the HI distri- 
bution is discussed. The main results are summarized in 
Sect.il 

2. SNe driven supershells: dynamics and 
morphology 

2.1. Model 

The equilibrium HI distribution is obtained by solving the 
steady-state momentum equation in cylindrical (r, z) co- 
ordinates in a fixed gravitational potential determined by 
the stellar and dark matter components. The self-gravity 
of gas is neglected. The dark matter density distribution 
(pd) is assumed to be that of a modified isothermal sphere 
(Binney & Tremaine ITM7)l 



Pd 



PdO 



(1 + r/r ) 2 



(1) 



where r is the radial distance from the galactic center, 
Pdo is the central density, and ro is the characteristic scale 
length. The stellar component is assumed to be vertically 
stratified (Spitzer EMU 

= p* sech 2 (z//i»), (2) 

where p, is the stellar density, p,o is the stellar midplane 
density, and is the vertical scale height of stellar distri- 
bution. 

The resulting equilibrium HI distribution is distinct 
for different dark matter and stellar distributions. Its ver- 
tical shape is neither Gaussian nor exponential (Cclnik ct 
al. I1979[) . however, it is better fitted by a Gaussian than 
by an exponent. We consider three different models with 
three different gas vertical scale heights. The parameters 
of the stellar and dark matter distributions for each model 
are shown in Table 1. Fitting the equilibrium vertical HI 
distribution of each model by a Gaussian, we obtain the 
1 a gas scale heights of /i=210, 370, and 500 pc at the 
galactocentric radius of 2 kpc. For comparison, the corre- 
sponding exponential gas scale heights are ft. C xp=230, 400, 
and 540 pc, respectively. 

Further, we assume that the size of SN-driven shells is 
much smaller than that of a galaxy. The radial dependence 



of the dark matter gravitational potential within the shell 
can therefore be neglected and the initial equilibrium con- 
figuration of gas becomes a function of only the distance 
above the midplane of a galaxy. In all models the initial 
gas velocity dispersion was taken as v g — 9 km s _1 . 

The energy of supernova explosions is released in a 
sphere with a radius of four zones. We use the constant 
wind approximation described in detail in Mac Low & 
Ferrara ( 1 9(22! • We convert the energy of each SN ex- 
plosion (10 51 ergs) totally into the thermal energy, be- 
cause in the present simulations we deal with large stel- 
lar clusters capable of producing a hundred supernovae. 
With such a number of SN explosions, the surrounding 
ISM will be quickly heated and diluted, making radiative 
cooling within the injection sphere ineffective. We choose 
the energy input phase to last for 30 Myr, which roughly 
corresponds to a difference in the lifetimes of the most and 
least massive stars capable of producing SNe in a cluster 
of simultaneously born stars. 

SNe generate a supersonically expanding wind that 
compresses the gas, thus creating a bubble filled with 
the hot ejected gas surrounded by a shell of compressed 
cold material. The gas dynamics is followed by solving the 
usual set of hydrodynamical equations in cylindrical coor- 
dinates using the method of finite-differences with a time- 
explicit, operator split solution procedure of the ZEUS- 
2D numerical hydrodynamics code described in detail in 
Stone & Norman (19521 ■ We have implemented the cool- 
ing curve of Bohringer & Hensler (1989) for a metallicity 
of one tenth of solar. We use an empirical heating func- 
tion tuned to balance the cooling in the background at- 
mosphere so that it maintains the gas in hydrostatic equi- 
librium and may be thought of as a crude model for the 
stellar energy input. Cooling and heating are treated nu- 
merically using Newton- Raphson iterations, supplemented 
by a bisection algorithm for occasional zones where the 
Newton-Raphson method does not converge. In order to 
monitor accuracy, the total change in the internal energy 
density in one time step is kept below 15%. If this condi- 
tion is not met, the time step is reduced and a solution is 
again sought. 

2.2. Qualitative description 

At the initial stages the remnant comes into a radiative 
phase relatively early, approximately at 

1/2 



10° 



yr, 



-^38 

n 

when its radius is only 



R~ 25 



L 



1/2 



pc, 



(3) 



(4) 



where L^s is the mechanical luminosity in 10 38 ergs s _1 re- 
leased by multiple SN explosions, n is the ambient volume 
density. This implies that the remnant is always radia- 
tive before it expands out to R ~ 2h and enters a break- 
through phase (i.e. the shell breaking out of the disk and 
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losing metal-enriched material to the intergalactic space). 
The breakthrough occurs at approximately 

t~l(T 3 (^j fc 4 yr, (5) 

where h is in parsecs. The further evolution shows a very 
fast drop of pressure inside the hot remnant and the 
"snowplow" expansion of a quasi-cylindrical shell of com- 
pressed cold material in the disk. Due to a sharp drop 
of pressure inside the remnant the expanding shell starts 
thickening and can no longer be treated as a thin shell. 
Figure la shows the temporal evolution of the inner and 
outer radii of the shell expanding in a face-on galaxy with 
a gas scale height of h = 500 pc and HI surface density of 
Ehi = 10.7 Mq pc -2 . The surface density of the swept-up 
gas is generally not constant across the shell, it reaches 
a maximum near the outer edge of the shell and falls off 
on both sides. We define the inner radius of the shell as a 
radial distance where Shi drops to its unperturbed value 
in the inner side of the shell. The outer radius is defined 
in the same manner. The mean radius is then obtained as 
an arithmetic average of the inner and outer radii. The 
mechanical luminosity of = 1 is assumed, which is 
equivalent to 100 SN explosions during the first 30 Myr. 
The cooling rate for T > 10 4 K is taken from Bohringer & 
Hcnsler ( 1989), for T < 10 4 K the cooling function rapidly 
declines to zero. It is seen that the shell remains thinner 
than a quarter or less of its mean radius only in the SN- 
driven expansion phase during the first 30-35 Myr. At this 
evolutionary stage, the mean radius of the shell is well de- 
scribed by the R oc (L/p) 1 / 5 i 3/5 law (Castor et al. H975) 
plotted by the dotted line in Fig. where p = 0.5 po 
and po is the midplane mass density of the ambient gas. 
The shell grows very fast after t = 35 Myr when a break- 
through is to occur. At the time of a breakthrough, the 
mean radius of the shell is R w 1.6 h. After the break- 
through the shell enters the "snowplow" phase and its 
mean radius grows approximately as R oc t 1 ^ (as shown 
in Fig. by the dotted line) . 

In case of a smaller gas scale height (h = 210 pc), the 
breakthrough occurs at t=13.5 Myr when R w 1.8 h. The 
mean radius of the shell is poorly described by both the 
spherical {R cx i 1/4 , McCray & Kafatos HUSTjl and cylin- 
drical (R cx i 1//3 ) thin-shell "snowplow" expansion laws 
as shown in Fig. ^ by the dashed and dotted lines, re- 
spectively. Thus, it is obvious that when an HI hole seen in 
face-on galaxies is fitted by the spherical thin-shell "snow- 
plow" expansion law R ~ i 1 / 4 , the resultant energy input 
into the hole from SN explosions may be largely overesti- 
mated. 

2.3. Dynamics versus vertical gas scale height 

At the early expansion stages, when the shells remain 
nearly spherical, their kinematical and dynamical char- 
acteristics scaled to the proper mechanical luminosity and 
local density are similar in different environments (see 
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Fig. 1. a) The solid lines delineate the inner and outer 
radii of the shell expanding in a galaxy with h = 500 
pc. The dotted line gives the radius of the shell as pre- 
dicted from the R cx (i/p) 1/5 £ 3/5 law for the SN-driven 
expansion at t < 35 Myr (Castor et al. I1975|l and from 
the R oc i 1 / 3 law for the "snow-plow" expansion at 
t > 35 Myr. b) The same for h = 210 pc. The dotted 
line gives the radius of the shell as predicted from the 
R oc i 1 ' 3 law for a cylindrical "snowplow" expansion at 
t > 13.5 Myr, the dashed line shows a spherical "snow- 
plow" R oc t 1 / 4 expansion. 



discussion in Castor et al. I1975f) . Therefore, one should 
explore the later expansion stages in order to determine 
the large scale properties of a parent galaxy. At the later 
stages, the expanding shells are radiative and their expan- 
sion velocities are only weakly supersonic so that the cor- 
responding Mach number is M £ 1 (in most face-on galax- 
ies the expanding HI shells are marginally supersonic, or 
even subsonic as in the case of Holmberg II, see discussion 
in Shchekinov et al. 2001 ). In such conditions the compres- 
sion factor of gas inside the shells (mass density contrast) 
is not high: for isothermal shock waves it is p/po = 7M 2 , 
where 7 is the ratio of specific heats, and for weakly su- 
personic shocks at the later stages (t ~ 30 — 40 Myr) it 
can be around 7. 
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A directly observed quantity is the HI column density. 
For an HI hole in a face-on galaxy created by consecutive 
SN explosions, the HI column density in a ring of com- 
pressed gas surrounding the hole is 



AT(HI) 



2y/2R 3 /'- 



(6) 



where R is the ring radius, AR is its thickness, no is the 
midplane density of ambient gas no = N (BI)/h, N (Rl) 
is the column density of the background (unperturbed) 
HI gas. Substituting AR — R po/3p and p/po ~ 7 into 
Eq. © , one obtains a linear dependence of the relative HI 
column density on its radius in a face-on ring 



/Chi = 



_ N(m) R 



iVo(HI) 



(7) 



with the slope inversely proportional to the HI vertical 
scale height of a parent galaxy. The highest relative col- 
umn density /Chi is reached when the breakthrough occurs 
at R ~ a h, where a is a slowly varying function of the 
vertical gas scale height. Simulations in Sect. 12.21 indicate 
that a w 1.6 for h = 500 pc, while for h — 210 pc a w 1.8. 
An expanding ring reaches maximal relative column den- 
sity /Chi ~ 4 only in galaxies with rather low vertical scale 
heights h < 200 pc, where the compressed shell always re- 
mains supersonic. After the breakthrough, the compressed 
gas in the shell starts to move into the cavity because of 
a sharp drop in pressure, and the column density in the 
ring decreases approximately as cx R . 

When the gas scale height is sufficiently large, the shell 
expansion becomes sonic before breaking out of the disk, 
the associated perturbation propagates with the sound 
speed and simultaneously the shell itself starts swelling 
into the cavity, which results in a decrease of the gas col- 
umn density. This occurs when 



R 



OOQ rl/2 ,1/2 y -l/2 

zou^ 38 n wo Zj 10 pc, 



"100 ^10 



(8) 



where h 



100 



h/100 pc, £10 is the unperturbed surface 
density of gas (S s o) i n units of 10 M Q pc~ 2 . For in- 
stance, when the gas scale height is h = 500 pc and the 
unperturbed surface density E s0 = 10-5 M Q pc~ 2 , the 
shell becomes sonic when its radius is R ~ 500 pc. As a 
consequence, the shell starts swelling into the cavity and 
the corresponding relative column density of the ring de- 
creases. 

Figure [3 shows the relative column density in the ring, 
/Chi, obtained in our numerical simulations for three differ- 
ent vertical gas scale heights and two values of the undis- 
turbed surface density S 9 o- The curves exhibit the char- 
acteristic behaviour described above, with the maximal 
values of /Chi being higher for the lower vertical gas scale 
heights. Specifically, for h — 210 pc the maximum value 
of /Chi is approximately 3.4 at the time when the shell 
starts swelling because of a sharp pressure drop after the 
breakthrough. For larger h the maximum value of /Chi is 
smaller. Note that /Chi is virtually independent of £ g o- For 
the total period of energy injection < 30 Myr and for the 



vertical gas scale heights < 600 pc, the maximum value of 
/Chi and the corresponding radius of a shell do not depend 
on the luminosity L, while in gaseous disks with h > 600 
pc the maximum value of /Chi and the corresponding ra- 
dius scale as L 1 / 2 . However, since a widely accepted value 
for the mechanical luminosity L — 10 38 ergs s _1 is typical 
for OB associations with a Salpeter IMF, one may expect 
that the curves depicted in Fig. [5] represent a universal 
relation and may be used to draw firm conclusions about 
the vertical scale heights of HI distributions in disk galax- 
ies. It is obviously seen that the /Chi versus R relation 
is characteristic for each vertical gas scale height only in 
the early phases when the shell has not yet broken out of 
the disk, while after the breakthrough /Chi can barely be 
distinguished among different h. Note in Fig. a range 
in which the relative column density in the ring varies in 
galaxies with low vertical gas scale heights. It is 1.4-3.4 for 
h = 210 pc, while /Chi is restricted to a narrower range of 
1.3-2 in galaxies with h — 700 pc. Thus, a narrow spread 
of observed column densities may indirectly indicate that 
the vertical gas scale height is large. Moreover, considering 
a universal and single- valued dependence of K-m(h, R) for 
R < 600 — 700 pc, the local gas scale height h can be in- 
ferred directly from the measurements of /Chi - Further, the 
measurements of /Chi for the individual HI holes located 
at different galactocentric distances allow for the determi- 
nation of the radial variations in the gas scale height. 

The increase in the mass of the dark matter gener- 
ally makes a gas disk thinner, i.e. the gas scale height h 
decreases as the dark matter mass increases. Assuming 
the local plane-parallel gas distribution, the total density 
Ptot(®T r ) m the disk at the galactocentric radius r and 
z = can be derived from h(r) = v g / AnGptot ( , r) 
(van der Kruit 198lJ. Provided that the radial distribution 
of the visible mass is known from independent measure- 
ments, one can obtain the radial distribution of the dark 
matter in the galactic plane. These are, of course, very 
approximate estimates; nevertheless, a similar approach 
applied to the dwarf irregular galaxy Holmberg I by Ott 
et al. H2001f) has yielded a dark matter mass within the HI 
content of the galaxy, the value of which is in agreement 
with that obtained in a more sophisticated numerical mod- 
eling by Vorobyov et al. (2004). Note, however, that the 
flaring and warping of the gas disk may further complicate 
this analysis. 

The velocity field of a shell expanding in a nearly face- 
on galaxy can additionally hint at the vertical structure of 
gas. At the initial stages far from breakthrough, the ver- 
tical component of gas velocity in an expanding shell, v z , 
scales as y/l — (r/i?) 2 , where r is the projected distance 
from the geometrical center of a shell and R is the shell 
radius. However, when the shell is near breakthrough, a 
substantial amount of gas constituting the shell becomes 
involved in strong vertical motion. In Fig. [3] we plot the 
mass-weighted v z as a function of the projected distance 
from the geometrical center of a shell for different vertical 
gas scale heights and different fixed shell radii. The solid 
lines represent the shells with a radius of R = 0.5 kpc, 
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Fig. 2. Relative column density in the ring versus its ra- 
dius: dashed, solid, and dotted lines represent h= 210, 500, 
and 700 pc, respectively; open circles correspond to T, g o = 
15.2 M Q pc" 2 , filled circles - £ 5 n = 10.7 M & pc" 2 . 



Fig. 3. Mass- weighted z-velocity profiles versus projected 
distance from the geometrical center of the shells. Solid 
and dotted lines represent the shells with the radii of 0.5 
and 0.7 kpc, respectively. Open squares - h = 210 pc, 
t = 18 Myr; open triangles - h — 370 pc, i=13.5 Myr; 
open circles - h = 500 pc, t = 13.5 Myr. Filled squares 
- h = 210 pc, t = 36 Myr; filled triangles - h = 370 pc, 
t = 31 Myr; filled circles - h = 500 pc, t = 29 Myr. The 
unperturbed HI surface density is T, g0 = 10.7 M Q pc~ 2 . 



while the dotted lines correspond to the shells with a ra- 
dius of R = 0.7 kpc. 

It is seen that the expanding gas is only weakly su- 
personic (the corresponding Mach number of 2-3) in the 
shells that are at the early phases of expansion, well before 
breakthrough (R = 0.5 kpc, h = 370 and 500 pc). On the 
other hand, the shell that approaches the breakthrough 
phase (R = 0.5 kpc and h — 210 pc) shows a steep veloc- 
ity increase up to v z = 60 km s -1 at the inner edge of the 
projected shell. The v z radial profiles remain qualitatively 
similar at different phases of the shell expansion, namely 
at R, = 0.5 kpc and R = 0.7 kpc. 

At the breakthrough phase, the radius of a shell is 
almost twice as large as the gas scale height h. From this 
point of view, one may expect that the most extended HI 
holes in face-on galaxies (such as HI holes number 2, 5, 8, 
10 and others in Holmberg II, Puche et al. 119921 ) whose 
z-velocity reveals such a behaviour at the inner edges can 
provide us with a direct measure of the local gas scale 
height. Specifically, scanning the HI radial velocity profiles 
within a hole with a beam smaller than its size (for the 
Ho II case, the beam must be less than 4" — 6") could 
help to identify nearly breaking-through shells, the radii 
of which could then serve as a direct measure of the local 
gas scale height. 

Out-of-plane SN explosions in nearly face-on galaxies 
can also provide quite a robust tool for the determination 
of the gas vertical scale height. In Fig. 0]we plot typical 
HI spectra of the shells produced by 100 successive SNe 
located at 100 pc above the midplane. The HI spectrum 



is obtained for a 1 kpc diameter beam centered on the ex- 
panding shell. Two model galaxies with different gas scale 
heights are considered: h = 500 pc and h — 210 pc. A 
one-dimensional velocity dispersion of er g = 3 km s -1 is 
assumed for the gas forming an expanding shell when con- 
structing the model HI spectra, while <j g = 9 km s _1 is 
adopted for the unperturbed gas. The HI spectrum of the 
shell expanding in a galaxy with h = 500 pc is character- 
ized by a typical two-humped structure. A smaller part 
of the accelerating shell expanding upwards (v z > 0) and 
a more massive part of the shell expanding downwards 
(v z < 0) with approximately constant velocity are clearly 
seen. The HI spectrum of the shell expanding in a galaxy 
with h = 210 pc has a single peak located at v z ^ 0. 
Velocity centroids of gas expanding downwards are obvi- 
ously smaller for lower gas scale heights h, because lower 
scale heights lead to a faster onset of breakthrough and, 
therefore, to a faster drop of pressure inside the SN-driven 
shell. 

In Fig. we plot the absolute values of the peak v z 
velocities of gas expanding downwards as a function of 
the corresponding HI intensities (see Fig. for three dif- 
ferent gas scale heights: h = 210, 370, and 500 pc. It is 
clearly seen that the smaller peak z- velocities are associ- 
ated with the shells expanding in galaxies with lower gas 
scale heights. A distinguishing feature of this dependence 
is that the peak expansion velocity is virtually indepen- 
dent of intensity in a relatively wide range covering the lat- 
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Fig. 4. The HI spectrum of an expanding shell produced 
by 100 successive SNe located 100 pc above the midplane 
in a galaxy with a) h = 500 pc, b) h = 210 pc. The 
HI spectra are numerically obtained for a 1 kpc diameter 
beam centered on the shell at three different times in Myr 
as indicated in each frame. 

est stages of expansion before and after the breakthrough, 
and varies approximately as 



peak 



0.9 



-0.24 



km ; 



(9) 



loopcy \ v ioopc / 

where z* is the location of SNe above the midplane. 

3. Effects of external ionizing radiation 

One of the most striking facts related to HI holes in dwarf 
face-on galaxies is the lack of an underlying stellar pop- 
ulation that could produce a strong energy input into 
the ISM. This tendency is obviously seen in Holmberg 
II (Ho II) where most extended HI holes are located out- 
side the stellar disk (Rhode et al. 119991 Stewart et al. 
2000| ■ Hence, the SN origin of those HI holes in Ho II is 
questionable. An alternative explanation of the HI hole 
formation connected with high-velocity HI clouds (HVCs) 
falling onto the disk seems to be excluded in the case of 
Ho II, because no HVCs are found nearby. In this case, 
Efremov et al. I|199 8) argue that a supernova in the merg- 
ing event of a compact binary system may be responsible 
for HI hole creation. However, a distinguishing feature of 
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Fig. 5. The peak z-velocity versus HI intensity for dif- 
ferent gas scale heights and locations of SN explosions 
above the midplane: dashed line - h = 210 pc, SNe at 
z* = 100 pc; dotted line - h = 370 pc, SNe at z* = 100 pc; 
solid line with filled squares - h = 500 pc, SNe at 
z* = 100 pc; solid line with open squares - h — 500 pc, 
SNe at z» = 50 pc. Different evolutionary stages of the 
shells are labeled in Myr. 



multiple HI holes on the periphery of Ho II is that they 
form quite a regular structure resembling a spiral wave 
(see Stewart et al. I2()00|l . and seem unlikely to be pro- 
duced by regularly placed supernovae explosions. In this 
connection it is worth mentioning the role of an exter- 
nal ionizing radiation field in shaping the HI distribution, 
particularly on the galactic periphery. 

The position of the boundary between HI and 
HII layers, Zi, in an exponential vertical gas distri- 
bution ionized by a given flux of external UV pho- 
tons / is determined by the following equation 2?tl = 
a r n,Q J exp(—2z/h)dz. The observed column density of 
atomic hydrogen, iV(HI) = 2uq J q ' exp(—z/h)dz, can then 
be found as 



JV(HI) - A^o(H) 



inlh 



(10) 



where a r is the hydrogen recombination rate and Nq (H) = 
2noh is the total column density of hydrogen (both neutral 
and ionized), Uq is the total midplane density of hydrogen. 
It is readily seen that for a typical value of the background 
ionizing flux of / ~ 10 6 photons cm~ 2 s _1 sr _1 and for 
h ~ 300 pc the column density of ionized hydrogen can 
reach 3 x 10 20 cm" 2 , nearly coincident with the observed 
in Ho II at the radii of ~ 6 arcmin (Puche et al. ^992 ) . 
This implies that the observed A?" (HI) may represent a 
small fraction of the total amount of hydrogen, and thus 
a relatively weak perturbation in no and/or in A^o(H) can 
cause a rather strong variation in the amount of neutral 
hydrogen. This effect increases for a flaring gaseous disk 



E. I. Vorobyov and Yu. A. Shchekinov: HI holes in galactic disks 



7 




R (arcmin) 



Fig. 6. Comparison of the total hydrogen column density 
and the observed HI in Ho II: solid line - observed HI 
(Puche et. al. 1992), dashed line shows total column den- 
sity JV (H) from (JTDJ with h = 625 pc, dotted line - iV (H) 
with h = 625 + 18.6(i? - l) 2 . 



with the scale height h growing outwards. In Fig. we 
plot the total column density of hydrogen, iVo(H), as a 
function of Ho II radius. A*o(H) is derived from the ob- 
served column density of neutral hydrogen TV (HI) for two 
models of the vertical gas distribution in Ho II: a plane- 
parallel distribution with h—625 pc and a flaring disk with 
h = 625 + 18.6(i?— I) 2 pc, where R is given in arcminutcs. 
In the outer regions at R £ 6 arcmin, the observed TV (HI) 
may represent a relatively small (less than 50 %) fraction 
of the total hydrogen content, and thus can be easily dis- 
turbed, ionized, and transformed into a hole of depleted 
HI column density. 

In principle, HI holes created by multiple SN ex- 
plosions can become radially asymmetric due to radial 
variations in the gas density, particularly at larger dis- 
tances from the galactic center where the flaring of the 
gas disk may dominate. For an exponential disk with 
the density distribution p oc exp(—z/h — R/Rq), the ex- 
pansion velocity of a shell at a given point varies as 
v s ex exp(z/2h + R/2Rq). In this case, the geometrical 
center of an HI shell can be shifted outward with respect 
to the position of SN explosions (Kalberla, private commu- 
nication). Typical scales in vertical and radial directions 
are of the order of hundreds and thousands of parsecs, re- 
spectively. For example, for the Ho II galaxy h ~ 600 pc, 
while the radial scale for the HI surface density i?o ~ 3 kpc 
(Puche et al. 1992, Bureau & Carignan, 2002). For face-on 
galaxies, direct measurements of the radial distributions 
of the volume density or flaring are not available, but one 
can assume that the radial scale lengths for them are of 
the same order as for the HI surface density. The corre- 
sponding deviations of the expansion velocities in radial 
5vr and vertical Sv z directions from a precise spherical 



shape can be estimated then as 5vr/5v z ~ h/Ro <C 1. In 
the case of Ho II, it becomes 5vr/5v z ~ 0.2. Hence, one 
may conclude that the expected radial asymmetry of a 
shell, and the corresponding displacement of the geomet- 
rical center of a shell with respect to the location of SN 
explosions, is less than 20% of the gas scale height. 

4. Summary 

Clustered supernova explosions in the disk of a face-on 
galaxy produce an expanding shell of compressed mate- 
rial, which is seen as an HI ring surrounding a central HI 
depression. In this paper we show that several physical 
characteristics of expanding shells are sensitive to the gas 
scale height and, hence, can be used for the determination 
of the dark matter content in face-on disk galaxies. 

• A functional dependence of the relative column den- 
sity of gas in the HI ring (/Chi) on the radius of the ring is 
found to be specific for a given gas scale height h, which 
makes it possible to infer h from observations of even a 
single HI ring. This is particularly true for the shells near 
breakthrough phase. 

• The vertical component of gas velocity in an expand- 
ing shell reveals a characteristic behaviour (particularly, at 
the stages close to a breakthrough), which can be used to 
trace the gas scale height h in galaxies with a sufficiently 
thin HI layer. 

• Out-of-plane SN explosions produce asymmetric ex- 
pansion, with a more massive part of the shell propagat- 
ing towards the denser regions of the disk. The velocity 
of this part of the shell is a single-valued function of the 
scale height h and the height above the midplane at which 
the SN explosions take place. 

• Some of HI holes, particularly in the outer regions 
of a galaxy, can be connected to the ionization of HI layer 
by extragalactic UV photons. 
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